Transition from Gaussian-orthogonal to Gaussian-unitary ensemble in a microwave 

billiard with threefold symmetry 
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Recently it has been shown that time-reversal invariant systems with discrete symmetries may 
display in certain irreducible subspaces spectral statistics corresponding to the Gaussian unitary 
ensemble (GUE) rather than to the expected orthogonal one (GOE) . A Kramers type degeneracy is 
predicted in such situations. We present results for a microwave billiard with a threefold rotational 
symmetry and with the option to display or break a reflection symmetry. This allows us to observe 
the change from GOE to GUE statistics for one subset of levels. Since it was not possible to separate 
the three subspectra reliably, the number variances for the superimposed spectra were studied. The 
experimental results are compared with a theoretical and numerical study considering the effects of 
level splitting and level loss. 

PACS numbers: 05.45.Mt, 11.30.Er 
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I. INTRODUCTION 

It is generally accepted that the eigenvalues of sys- 
tems, whose classical counterparts are fully chaotic, be- 
have similarly to those of random matrix ensembles ||, |) . 
This conjecture also has some theoretical foundations 
||, ^, [^]. If we disregard spin there are two random ma- 
trix ensembles, with two different behaviors under time 
reversal. Systems with time reversal invariance (TRI) 
are described by an ensemble of real symmetric matri- 
ces, corresponding to the invariance of the Hamiltonian 
under complex conjugation. This ensemble is known as 
the Gaussian Orthogonal Ensemble (GOE) and should 
be contrasted with the ensemble of all Hcrmitian matri- 
ces, known as the Gaussian Unitary Ensemble (GUE), 
which usually applies when TRI is broken. 

An additional issue arises if discrete symmetries are 
present. In this case, it has been assumed that eigen- 
values belonging to different irreducible representations 
of the symmetry group are statistically independent and 
obey GOE or GUE spectral statistics, depending on 
whether TRI holds or not. These assumptions have re- 
ceived considerable support both from theoretical con- 
siderations and from numerical simulations. However, in 
recent papers || [7| , it was suggested that an anomalous 
situation arises when a time-reversal invariant system has 
a discrete symmetry, such that TRI is broken within one 
of the irreducible subspaces of the point group. This is 
equivalent to the fact that the symmetry group has non- 
selfconjugate representations. In this case, TRI induces a 
degeneracy between different (conjugate) irreducible rep- 
resentations of the symmetry group (Kramers degener- 
acy). It was suggested that the eigenvalues belonging 
to these degenerate subspaces should have the statistical 
properties of the GUE rather than the GOE. This might 
be expected from the fact that complex conjugation actu- 
ally maps one irreducible representation onto the other, 



rather than leaving it invariant. In this sense, one may 
indeed argue that TRI is violated within the subspace 
belonging to either representation, and is only restored 
through the combined presence of both representations. 
A more complete argument is presented in [Q. A semi- 
classical view of this question is given in || . At this point 
we should add a note of caution. This entire argument 
applies to the largest symmetry group of the system, and 
relies on the fact that this group is known. Deviations 
from the expectations would indicate the presence of a 
larger unknown symmetry group. On the other hand, 
symmetry breaking, which is usually present in experi- 
ment, will also affect the results. In the present paper we 
will address both points in a simplest example. 

Let us now consider a system with threefold symmetry 
but with no twofold symmetry axis, i. e. C3 symmetry 
This means that the wave functions in polar coordinates 
can be put in the form 



ipn,a(r,' 



J2 <fe,nMe 4(3fc+CT)0 



(1) 



where n runs over all the integers and a takes the values 
and ±1. Because of TRI, it is clear that the eigenval- 
ues corresponding to ij) n ,i and VV-i are degenerate, and 
we expect them to show GUE statistics. If in addition 
there is a twofold symmetry axis leading to C% v sym- 
metry, the two conjugate representations combine to a 
two-dimensional irreducible representation which is self- 
conjugate. Thus we expect GOE statistics as well as 
degeneracy. 

There is one previous microwave experiment on a bil- 
liard with C3 symmetry performed by the Richter group 
H . In the data analysis the authors faced the problem to 
separate the subspectra belonging to different irreducible 
representations. Since it is impossible to realize a C3 
symmetry exactly in experiment, the doubly degenerate 
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eigenvalues belonging to the a = ±1 representations are 
split into doublets. In Rcf. |IJ this doublet splitting was 
used as a tool to separate the subspectra by classifying 
each pair of eigenvalues with a distance smaller than some 
given threshold value as a split doublet. The technique 
worked quite well, but there was a misidcntification of 
about 5 percent of the levels. Nevertheless, the authors 
were able to see the GOE and GUE signatures in the 
respective subspectra. To improve the agreement with 
random matrix predictions, the authors calculated the 
spectrum numerically and used the results of the calcu- 
lation to relabel the misidentified eigenvalues. 

In the present work a billiard with threefold symmetry 
is studied, which has the additional feature to allow a 
change of symmetry from Cs, v to C3 by turning an in- 
sert placed in the center of the billiard, thus allowing to 
study the transition from GOE to GUE behavior of the 
(7 = ±1 subspectra. In view of the problems to separate 
the respective subspectra satisfactorily, we refrained from 
any attempt in this direction. Instead we concentrated 
on the spectral statistics for the superimposed spectrum 
and compared the results with numerical and analytic 
calculations based on random matrix theory (RMT). 

A formula for the number variance is derived assum- 
ing that the degeneracy is lifted by a random off-diagonal 
matrix element between the two states of the doublets. 
It describes fairly well the long range properties of the 
spectra completed by means of level dynamics. The spec- 
tral short range behavior is sensitive to the details of the 
completion procedure while the long range behavior is 
not. This is why the difference between the GUE and 
the GOE case can still be detected very reliably in the 
long range regime. 

If we consider the individual spectra there is a loss 
of roughly 20% of the levels. The number variance of 
the uncompleted spectra were compared with numerical 
RMT simulations. Several ways to simulate the loss of 
levels have been applied, and in all cases the two sym- 
metries remained distinguishable at a loss rate of 20%. 

In the following section the experimental setup and 
the details of the data analysis are presented. Further 
the gradual transition from the C% v to the C3 symmetry 
will be observed in the number variance of the experi- 
mental spectra. In section III the results for the com- 
pleted spectra will be compared with calculations taking 
into account the break of the C3 symmetry. Finally, in 
section IV we will discuss the results for individual spec- 
tra with missing levels; these will be complemented with 
numerical RMT simulations. 



II. EXPERIMENTAL TECHNIQUE 

The main features of the experimental technique have 
been presented elsewhere |l(J. Therefore we concentrate 
on the aspects of relevance in the present context. Figure 
[l] shows a sketch of the billiard used in the experiment. 
The interior insert can be rotated versus the outer part 




FIG. 1: Sketch of the microwave billiard used in the exper- 
iment at angle <f> — 0°. The positions of the antennas are 
marked by crosses. 

thus allowing a spectral level dynamics measurement in 
dependence of the rotation angle 4>. To avoid a break 
of symmetry due to the presence of the antenna, three 
antennas were placed at symmetry equivalent positions, 
two of them being terminated by 50 f2 loads. Spectra 
were taken for angles from <f> = 0°, corresponding to the 
geometry with C$ v symmetry shown in the figure, to <f> = 
20° in steps of one degree. All boundaries, both of the 
insert and the outer part, were curved, and the corners 
of the outer boundary were removed to avoid families 
of marginally stable orbits such as bouncing balls. This 
was checked by means of a Poincare plot of the classical 
motion of a particle in the billiard. 

The height of the billiard was h = 10 mm, i. e. the 
system was quasi-two-dimensional for frequencies below 
f max = c/2h = 15 GHz. Microwave reflection spectra 
were taken with a Wiltron 360B vector network analyzer 
in the frequency range 0.5 to 7 GHz with a resolution 
of 0.5 MHz. In the studied spectral range one expects 
to find altogether 400 eigenvalues for each angle <f>. This 
number is obtained from the Weyl formula with area and 
circumference term 

n(k) = 4-k 2 -^-k. (2) 

n(k) is the mean integrated density of states, where 
A « 0.258 m 2 and S « 3.59 m are area and circumfer- 
ence of the billiard, and k — 2nv/c is the wave number. 
The curvature of the billiard yields a correction term of 
the order of 1. The number of states identified for a 
fixed angle falls approximately 20% short of the Weyl es- 
timate. It is possible to recover all missing eigenvalues by 
studying the dynamics of the spectra as a function of the 
angle. This is demonstrated in Fig. || where the positions 
of the missing resonances are marked by diamond-shaped 
symbols. 

We begin with a presentation of the number variance 
in dependence of the insert angle, both for the completed 
and uncompleted spectra. In both cases the spectra were 
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FIG. 2: Part of the spectral level dynamics for the billiard 
with threefold symmetry with the orientation angle <f> of the 
central insert as the level dynamics parameter. Angle (j> = 0° 
corresponds to C*3„ symmetry. Only the resonances marked 
by cross-shaped symbols were really found in the experiment. 
At the positions marked by diamond-shaped symbols missing 
resonances were supplemented. 

unfolded to a constant mean density of states using a sec- 
ond order polynomial fit to the experimental integrated 
density of states. The number variance E 2 (L) has been 
obtained by moving the interval of length L through each 
spectrum on a fine grid to not loose any information. 
Consequently, the smoothness of S 2 (L) does not reflect 
its statistical quality, since the intervals are not statisti- 
cally independent. 

Figure || shows the number variance E 2 (£) for differ- 
ent angles </>, where the results for <fi = 0°, 1°, <fi — 2°, 3°, 
cf> = 4°, 5° and <p — 6°, 7° have been averaged. We see a 
gradual transition from the C% v to the C3 case with in- 
creasing <f). The line with long dashes shows an average 
over <j> — 8° . . . 19°, corresponding to the C3 case. The 
GOE-GUE transition is clearly seen, both in the com- 
pleted and the uncompleted spectra. In the following we 
will only discuss the two extreme cases, where for the C% v 
case the results for <\> = 0°, 1° will be averaged and for 
the C3 case the ones for <f> = 8° . . . 19°. 

III. RESULTS FOR THE COMPLETED 
SPECTRA AND THE DOUBLET SPLITTING 

We recall that for <fi = 0° we expect a superposition 
of one GOE spectrum and another doubly degenerate 
GOE spectrum for the ideal system, whereas for angles 
4> sufficiently different from zero a superposition of one 
GOE spectrum with a doubly degenerate GUE spectrum 
should be observed. The short range behavior is domi- 
nated by doublet splitting, which we have no reliable ex- 
perimental handle on. We therefore concentrate on the 
long range behavior. In particular we consider the long 
range part of the two-point correlation function. The 
number variance is ideally suited for this purpose because 
of the clear signature it shows for the difference between 
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FIG. 3: Number variance S 2 (L) for different angles <j> for 
the completed (a) and the uncompleted (b) spectra. The 
results for <j> = 0°, 1° (upper solid line), <j> = 2°, 3° (dashed), 
4> = 4°, 5° (dash-dotted), <f> — 6°, 7° (dotted), were averaged. 
The lower solid line shows an average over = 8°... 19°. 

the GOE and the GUE case. 

For the superposition of M strictly non-interacting, 
equally weighted subspectra the number variance can be 
written as 

M 

^total(^) = (if-k) ■ 

m— 1 

In the case of a threefold symmetry, the three subspectra 
are not independent; instead two of them are degenerate 
and thus the covariance of this doublet spectrum leads to 
a factor 2 in the number variance of their superposition: 

£ t 2 otal (L) = £ 2 (iL)+4£ 2 (iL) . (3) 

Here E 2 is the number variance of the singlet spectrum 
being a GOE spectrum in both cases, whereas E| is the 
number variance of the doublet spectrum which may be 
a GOE or a GUE spectrum, depending on the symmetry 
of the system. 

In the analysis of the completed spectra, the first 70 
eigenvalues of each spectrum have been omitted to avoid 
non-generic features in the statistics. The results for the 
number variance of the experimental spectra have been 
averaged over the rotation angles = 0, 1° in the C% v 



4 



3.0 




■ I ' I ■ I ■ I ■ 1— h 




FIG. 4: Number variance S 2 (L) of the completed spectra for 
the mirror-symmetric (a) and the non-mirror-symmetric (b) 
case. The experimental results (solid line) are compared with 
the theoretical prediction for the ideal system (dotted line), 
and with an analytical result, where the doublet splitting was 
taken into account (dash-dotted line). 



case and <f> = 8.. 19° in the C3 case. This implies that the 
statistical fluctuations of the number variance are much 
larger in the former case. 

Figure |] shows experimental data with solid lines for 
the Csv (a) and the C3 (b) case, respectively, while the 
dotted lines give the ideal result of Eq. |[ It is not sur- 
prising that the agreement is poor, since we already no- 
ticed in the level dynamics that the doublet spectrum 
was rather strongly split, indicating a sizeable break of 
the C3 symmetry of the billiard. On the other hand, 
we clearly see that the two number variances differ, and 
that the effect is quite large. This was the motivation 
for further theoretical considerations of the symmetry 
breaking. A detailed discussion of this issue follows in 
the appendix, where an expression for the corresponding 
number variance is derived. Here only the main points 
shall be outlined. 

The case of the C3 symmetry is the more delicate one. 
The basis functions spanning the irreducible subspaces 
are complex, but due to the degeneracy we can always 
choose real functions. The appendix gives formal expla- 
nations, but here it may suffice to say that the actual 
measurements are on real fields and therefore it is clear 



that the smallest perturbation that does lift the degener- 
acy must cause real linear combinations. Therefore it is 
not surprising that we have to use degenerate perturba- 
tion theory for a real symmetric matrix. Furthermore it is 
quite clear that the perturbation matrix elements should 
be Gaussian distributed. While the mechanical imper- 
fections causing the breakdown of the symmetry are the 
same for all functions, the unperturbed eigenfunctions 
involved are those of a chaotic system and thus in good 
approximation random functions that should yield Gaus- 
sian distributed matrix elements. Based on these simple 
considerations it is clear that the level splitting is Gaus- 
sian distributed, too, and linear in the perturbation. In 
the appendix we obtain an expression for the effect of 
such a level splitting on the long range part of the num- 
ber variance, and we also show that effects of three level 
terms are small and thus effects of the symmetry break- 
ing on the number variance for larger distances are of 
second order in the perturbation. 

This result is of central importance since it explains 
why the difference in the number variances persists in 
spite of the fairly large doublet splitting. A parameter 
A = 0.125 for the width of the distribution of the pertur- 
bation matrix element and thus for the doublet splitting 
(see Eq. A5 ) , is consistent with the splitting found in the 
experiment. It corresponds to a mean level splitting of 
about 0.85 in units of the mean level spacing. If we now 
use the results of the appendix for the number variance, 
Eq.|Al"(i|, both for the GOE and the GUE case, we find 
good agreement with the experimental results, as can be 
inferred from the dash-dotted line in Fig.Q. The small 
deviations are consistent with the statistical uncertain- 
ties seen in skewness and excess (not shown). We there- 
fore have a clear understanding of the role played by the 
splitting of the degeneracy and why it does not affect the 
difference between the two cases even if the symmetries 
are only approximate. 



IV. MISSING LEVELS AND INDIVIDUAL 
SPECTRA 

Completing the spectra by means of spectral level dy- 
namics we were able to come to a quantitative theoretical 
understanding of the GOE-GUE transition observed in 
the experiment. The completion of spectra is definitely 
a legitimate approach, but in general only isolated spec- 
tra for fixed systems are available. Therefore we felt the 
need to check the robustness of the signatures of the tran- 
sition by considering what we can learn from incomplete 
spectra with as many as 20% of the levels missing. 

Numerical simulations have been performed as follows. 
We produced ensembles of 1000 spectra of dimension 500 
to create the GOE or the GUE doublets for the C% v 
and C3 case respectively. Then we lifted the degener- 
acy randomly using a Gaussian distribution with a width 
of about 12.5% of the average level distance of a single 
spectrum. These spectra have been superposed with an 
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independent GOE spectrum of the same strength. 

In the next step we tried to simulate the experimen- 
tal loss. There are two sources of missing levels. First, 
close pairs can no longer be resolved if their separation 
is smaller than the line width. Second, a level will be 
missed if incidentally a node line of the associated wave 
function will be just at the position of the antenna. This 
mechanism leads to a loss that is uniformly distributed 
along the energy axis. 

So we skipped randomly 7% of the eigenvalues to sim- 
ulate the global loss due to the position of the antenna. 
Then we compared the distance between two adjacent 
levels with a Gaussian distributed random number whose 
width is chosen such that around 13% of the doublets are 
destroyed. From these spectra only 30% of the levels at 
the center of each spectrum are taken into account for 
the further analysis in order to avoid edge effects. The 
ratio of uncorrelated and correlated loss has been cho- 
sen to get a fair correspondence between experiment and 
theory for S 2 at small distances and is consistent with 
values found in previous works (see e. g. Ref. |Tl[| ). 

The numerical as well as the experimental results for 
the nearest neighbor distribution for the case of C3 sym- 
metry are presented in Fig.|^. A good qualitative agree- 
ment between the two curves is found, although the first 
bins are underestimated in the numerical data resulting 
in an overshoot of its maximum. In Figure ^ the exper- 
imental results for both symmetries are compared with 
the numerics. In both cases we get a rough agreement up 
to L = 3. For larger distances the slope of the numerical 
curve is by far too large which indicates that we have 
somehow destroyed long range correlations still present 
in the experimental data. On the other hand, the relia- 
bility of the experimental data for large L is questionable 
because of the missing levels. Therefore it seems unrea- 
sonable to fit parameters like the percentages of the two 
sources of level missing. In any case, both the experimen- 
tal results as well as the numerical calculations permit a 
clear distinction between the two symmetry cases. 



CONCLUSIONS 




FIG. 5: Nearest neighbor distribution for the case of C3 sym- 
metry. The black line shows the results for the experimental 
data with 20% loss, the grey line shows the numerical simu- 
lations. 




FIG. 6: Number variance S 2 (L) of the experimental data 
with 20% loss (solid line: Cz v , dash-dotted line: C3 case) in 
comparison with the numerical simulations (dashed line: C3,,, 
dotted line: C3 case). 



We have presented experimental results for a quasi- 
two-dimensional microwave cavity with a movable insert 
permitting to change the discrete symmetry of the cavity 
from C% v to C3 . According to theory, this change should 
modify the spectral statistics in the ideal case from a su- 
perposition of two GOE spectra, one of them degenerate, 
to a superposition of a GOE spectrum with a degener- 
ate GUE spectrum. We find that the level splitting due 
to imperfections in the C3 symmetry amounts to about 
0.85 of the mean level spacing, and we face up to 20% 
loss of levels in an isolated spectrum. From spectral level 
dynamics by rotating the cavity's insert we can recover 
the entire spectra, with some errors in the short range 
behavior of correlation functions. Nevertheless, the long 
range part of the number variance still clearly displays 



the difference between the two cases. From a theoretical 
analysis of the symmetry breaking in terms of degenerate 
perturbation theory, this behavior becomes comprehen- 
sible, since the level splitting is linear in the perturbation 
while long-range effects are quadratic. 

Even more importantly, we have shown that the spec- 
tral statistics of spectra with large loss of levels can still 
provide usable information. By postulating two types of 
losses, a statistical one due to the fact that the antenna 
may be close to a node of the wavefunction and a sys- 
tematic one due to overlap of close-lying levels, we can 
explain the experimental results qualitatively. 
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APPENDIX A: NUMBER VARIANCE FOR 
PERTURBED C* 3 SYMMETRY 



perturbations as long as the Vi are small with respect 
to the third eigenvalue. If on the other hand the third 
eigenvalue is so close to the degenerate doublet that the 
V. are of the order of one, the situation becomes more 
complex and Cardan's formula must be used. However, 
the splitting will still be linear in the Vi, so that the 
qualitative behavior should not be greatly affected. The 
probability for one doublet to come close to another one 
is very small, since they experience GUE level repulsion. 

We now turn to the calculation of the spectral form 
factor, and of the number variance for a weakly split GUE 
doublet spectrum. Let us take the e k to be an arbitrary 
reference spectrum with mean level spacing equal to one. 
The split doublet spectrum is generated by 



Ek,<r = 2e k + aa k 



(A4) 



In this appendix we present a detailed treatment of 
the effect of symmetry breaking. We first need to spec- 
ify a reasonable RMT model for this system. To this 
end let us consider an unperturbed Hamiltonian matrix 
Hq representing the billiard with perfect C 3 symmetry. 
The Hilbert space can be separated into three irreducible 
parts, one invariant under the group action, the other 
two spanned by clockwise and counterclockwise travelling 
waves. Because of TRI the eigenfunctions in the latter 
two subspaces have pairwise degenerate eigenvalues, i. e. 
they display a Kramers degeneracy. If we wish to display 
TRI rather than the C3 invariance we can alternatively 
use a real basis for the doublet space. Any perturbation 
due to imperfections of the billiard can only break the C3 
symmetry but not TRI. As we wish to introduce an RMT 
model for the perturbation that conserves TRI, it is ob- 
vious that this can be implemented by a real symmetric 
matrix in the real basis. 

From this picture the effect of the perturbation is im- 
mediately clear: Since we have degenerate eigenvalues, 
we must look at these separately. Then we may ask what 
happens if there is a third eigenvalue in the vicinity. We 
therefore consider the matrix 



.4 



V 1 V 2 ' 
Vi V 3 

v 2 v 3 1 , 



(Al) 



where we have scaled and shifted the two degenerate 
eigenvalues to zero, and the perturbing third eigenvalue 
one. The eigenvalues A^ obey the relations 

Ai + A 2 + A 3 = 1 

A 2 + A 2 + A 2 = 1 + 2(V? + Vi + V 3 2 ) (A2) 
AiA 2 A 3 = -V? + 2ViViVa. 
These equations are satisfied to second order in the Vi by 

Ai = Vx 

A 2 = -Vl (A3) 
A 3 = 1. 

The shift in the two degenerate eigenvalues, being linear 
in the off-diagonal element V\ , thus dominates all other 



where a runs over ±1 and the factor 2 retains the mean 
level spacing of one. The a k are taken to be independent 
Gaussian random variables distributed according to 

p(flH \/^ eXP ^ a2/(2A) )- (A5) 



We are now going to compute the form factor 

N 

2X V 



*(*) = 5*F E £ eM*-.--*U0« (A6) 



a,a' = ±l k.l = l 



from the corresponding quantity ko(t) for the e k . Sep- 
arating the right hand side into a part for k = I and 
another for k ^ I and averaging over the a k yields 



k(t) = l + er m At -2e~ 47r At +2/c (2t)e" 47r ^ . (A7) 

The number variance S 2 (L) can be expressed as (see e. g. 
Ref. ffj) 



£ 2 (L) = — / dtk{t) 
Jo 



sin irLt 
t 



(A8) 



Thus, from Eq. A8 we get the S 2 (L) for the doublet part 



of the spectrum by inserting the GUE form factor 

k (t) = \t\ (\t\ < 1) 
= 1 (1*1 > 1) 



(A9) 



into (A7). From this one finally obtains for the number 
variance of the superimposed spectrum 



£ 2 otQi (L) = £ 2 (§L)+£ 2 , 



(A10) 



where S 2 (L) is the number variance of the GOE singlet 
spectrum. 

The C 3v case yields a similar result if we substitute the 
GOE form factor rather than the GUE form factor for k 



in Eq. A7 
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